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Abstract. Wc consider constrained Hamiltonian systems in the framework of Dirac's 
theory. We show that the Jacobi identity results from imposing that the constraints are 
Casimir invariants, regardless of the fact that the matrix of Poisson brackets between 
constraints is invertible or not. We point out that the proof we provide ensures the 
validity of the Jacobi identity everywhere in phase space, and not just on the surface 
defined by the constraints. Two examples are considered: A finite dimensional system 
with an odd number of constraints, and the Vlasov-Poisson reduction from Vlasov- 
Maxwell equations. 
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1. Introduction 

Imposing constraints on a Hamiltonian system is routinely done using Dirac's theory [Tj 
[21 El S] - The key point is to compute the matrix C of the Poisson brackets between 
two constraints (i.e., whose elements are C nm = $ m } where $ n (z) = are the 
constraints) and to invert this matrix. Under this hypothesis, it has been shown that 
the following bracket, called Dirac bracket, 

{F, G}* = {F, G} - {F, <Z> n }D nm {$ m , G}, (1) 

where D = C _1 , is a Poisson bracket. The technical difficulty is to prove the Jacobi 
identity, and this has been done by Dirac in Ref. [TJ, and his proof relies heavily on the 
invertibility of C (see Eq. (59) which results from Eq. (35) in Ref. pQ, or Appendix B 
in Ref. [5]). As a consequence of the fact that D is the inverse of C, the constraints are 
Casimir invariants, i.e., {$„, G}* = for any observable G. 

What if the matrix C is not invertible? This could happen for instance if there is an 
odd number of constraints (since C is antisymmetric) or if the constraints are reducible or 
redundant [HI [TJ [8] . The purpose of the present article is to address this question. Our 
main result is to show that the Jacobi identity is a property which results from imposing 
that the constraints are Casimir invariants, regardless of the invertibility of C. This 
result holds for canonical or non-canonical Hamiltonian systems (see Refs. [TUl [TTJ [12] 
for an introduction to non-canonical Hamiltonian systems). We notice that the proof 
we provide shows the validity of the Jacobi identity everywhere in phase space, and not 
just on the surface defined by the constraints as found in Refs. [6j [TJ [8]. 

We consider a finite dimensional Hamiltonian system whose variables are denoted 
z = (zi, z 2 , ■ . . , zn). It is given by a Hamiltonian H(z), a scalar function of the variables, 
and a Poisson bracket written as 

where the Poisson matrix j(z) is such that the bracket ([2]) is antisymmetric and satisfies 
the Jacobi identity [in addition to the bilinearity and the Leibnitz rule which are already 
ensured by the form of the bracket (121)]. 

We impose a set of M < N — 2 constraints $ m (z) = for m = 1, . . . , M which 
are scalar functions of the variables z. We consider brackets of the form (TTJ with an 
antisymmetric matrix D which is not necessarily the inverse of C (and consequently 
these brackets are not not Poisson brackets in general). The matrix associated with the 
bracket (TT]) is given by [13] 

J* = J-jQ f DQj, (3) 

where the matrix Q has elements Qm = dQ n /dzi. The matrix D is chosen such that the 
constraints are Casimir invariants. This leads to the following condition on D : 

jQ t (l-DC) = 0, (4) 

where 1 is the M x M identity matrix and C = QjQ} . A first situation is when C is 
invertible, and hence a possible solution to Eq. fll]) is D = C _1 . This is the main case 
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considered in the literature, for which the Jacobi identity is proved in order to ensure 
that the Dirac bracket is a Poisson bracket [TJ [5] . 

However in many instances, the matrix C is not invertible, and, to the best of our 
knowledge, there is no proof of the Jacobi identity in these cases which holds everywhere 
in phase space. There are two properties which are strongly dependent on the particular 
choice of matrix D, which are (i) the Jacobi identity and (ii) the fact that the constraints 
are Casimir invariants. The purpose of this article is to show that (ii) implies (i), or in 
other words that imposing that the constraints are Casimir invariants, i.e., Eq. (jl]), is 
sufficient to ensure that the Dirac bracket ([1]) is a Poisson bracket. The proof is done 
in Sec. El In Sec. [3] we illustrate the computation of the Dirac bracket in cases where 
C is not invertible, with two examples: a finite-dimensional Hamiltonian system with 
an odd number of constraints, and the Vlasov-Poisson reduction from Vlasov-Maxwell 
equations. 

Before going into the proof of the Jacobi identity, we would like to address two 
questions: 

1) Does the Jacobi identity for the bracket ([I]) imply that the constraints are Casimir 
invariants, i.e., does (i) implies (ii)7 The answer is no, and we provide a counter example 
below. Consider a Poisson matrix J written in block form as 



J 



where C and J satisfy individually the Jacobi identity so that the matrix J too. The 
bracket (CQ) is characterized by a matrix 



J, 



for Q = (l, 0). Assuming that C is invertible, we choose D = C _1 (l — A) so that 
C = C(l — DC) = AC which satisfies the Jacobi identity inherited from C. In order to 
have the constraints as Casimir invariants, C must vanish. Therefore the bracket (flj) 
satisfies the Jacobi identity in this case, but does not have the constraints (zi, . . . , Zm) 
(where M is the number of columns of C) as Casimir invariants. 

2) Given that Eq. (jl]) might have more than one solution, does it lead to different 
expressions for the Dirac bracket? The answer is no. We consider D a solution of Eq. (j3J). 
We notice that any matrix D = D + with JQ) C = also satisfies Eq. (j4j). The Dirac 
bracket is obtained from the Dirac projector [13J V — 1 — JQ^DQ, i.e., J* = VJV^. If 
we consider the other projector associated with D, i.e., V = V — JQ) Q, then we show 
that 

pjpt = -pj-p\ 

where we use the identity VJQ} = 0. This identity ensures that the Dirac bracket is 
unique, even if there might be more than one solution to Eq. pi) . 
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2. Proof of Jacobi identity 

A proof of a weak version of the Jacobi identity, i.e., the validity of the Jacobi identity 
on the surface defined by the constraints, has been detailed, e.g., in Refs. [6J [TJ, [8] . It is 
based on showing that 

{F, {G, #}*}* « {F', {G', if'}}, 
where F' = F — {F, § n }D nm Q m in order to deduce that 

{F, {G, H}*}* + {H, {F, G}*}* + {G, {H, F}*}* « 0. 

Here we show that the Jacobi identity holds everywhere in phase space, i.e., the weak 
equality can be made a strong one, i.e., 

{F, {G, H}*}. + {H, {F, G}*}* + {G, {H, F}*}, = 0. 

As a consequence, even in the case where the matrix C is not invertible, the Dirac 
bracket ([I]) [if it can be constructed using Eq. (J3J)] is a Poisson bracket. 

First we perform a local change of coordinates such that the new variables are the 
constraint functions. This can be done at least locally under the assumption of the 
change of coordinates. In other terms, we assume that $fc(z) = z k for k G [1,M]. We 
assume that J satisfies the Jacobi identity. In the variables (<&, w) the Poisson matrix 
is expressed by blocks 

--(:?)■ 

The Poisson matrix associated with the bracket (CEJ) is given by 

/ C(l - DC) -(1 - CD)B t \ 
I H(1 :;C! J + BDB f J' 

We assume that D is chosen such that the constraints are Casimir invariants of J*. This 



condition becomes: 

C(1-DC) = 0, (5) 
B(1-DC) = 0. (6) 

The Jacobi identity for J translates into four sets of equations 

CuViC jk - B H d l C jh + O m = for k e [1, M], (7) 



CjiViBki — CuViBkj + BudiBkj — B\jd\Bki + BkiViCij 

+ JkidiCij = for i,j e [1, M] and k e [1,N — M], (8) 
BjiViB ki — B k iViBji + JjidiB ki — J k [diBji + CuViJj k 

-B H diJ jk = for i E [1, M\ and j, k <E [1,N — M], (9) 
Bu^iJ jk + JudiJ jk + {ijk) = for k G [1, N - M), (10) 



Casimir invariants and the Jacobi identity in Dirac's theory 



5 



where Oajk) designates the terms obtained by circular permutations of the indices 
(i,j,k), V/ = d/d^i and di = d/dwi. So when the index I is involved with V, the 
implicit sum runs from Z = 1 to M. When it is involved with d, the implicit sum runs 
from 1 = 1 to N — M. 

Given Eqs. (JHJ)-(|EJ), the Jacobi identity for J* reduces to 

(J + BDB + )«^(J + BDB*) jk + 0(y fe) = 0. (11) 

The aim is to use Eqs. (!7|)- (|T0|) together with Eq. (Q in order to prove Eq. ({TTI) . 
Equation (TTTT) can be decomposed into three sets of terms : 

S ijk = JiAJjk + {BDB%d l J jk + 0( ijk ), (12) 
T ijk = J il d l (BDB^) jk + {ijk ), (13) 
U ijk = (BDB t ) li ^(BDB t ) ifc + O m . (14) 

Here we notice that all indices i, j, k belong to [1, N — M]. Using Eq. ( flQj) . the 5" terms 
can be rewritten as 

Sijk = —BaViJjk + B im D mn Bi n diJjk+ 0(ijk) ■ 

By rewriting B[ n diJjk using Eq. fl9]), a cancellation is obtained from Eq. (Q, and the S 
terms are rewritten as 

Sijk = B im D im {BjiV iB kn — B k iViBj n + JjidiB kn — J k idiBj n )+ 0(ij k ) ■ 

By noticing a cancellation in the terms BDJdB in S and T (using a circular permutation 
of the indices (i, j, k) and the antisymmetry of D), we obtain 

Sijk+Tijk = B im D mn (Bji'ViBkn—Bki'ViBj n )+JiiBj m BkndiD mn + 0(ijk) -(15) 

Concerning the U terms, we decompose them into two parts : 

^ ijk Bi m D mn Bi n DpqBjpdiB k q -\- Bi m D mn Bi n DpqB k qOiBjp-\- (3(ijk): 

L \jk Bi rn D mn Bi n BjpBkqdiDpq-\- O(ijfc) • 

The second term of is rewritten as BjpDp q D mn B im Bi q diB kn using a circular 
permutation of the indices (i, j, k). Therefore, we have 

Uijl = Bi m D mn Dp q Bjp(Bi n diB kq — Bi q diB kn )+ 0(ij k ), 

where we have used the antisymmetry of D. From Eq. ([8]), is rewritten as 

Uijk = ~ BimD mn (BT}C)ji¥lBkn + Bj m D mn (BDC)uVlBkn 

Bi m BjpBklD mn Dpq^ iCnq B i m B jpB mn B) p q J k [d\C nq -\- 0(jj'fc)(i6) 

From Eq. ([6]) and using a circular permutation on the indices (i, j, k), the first line of the 
previous equation cancels with the BDBVB terms in Eq. ( fl5l) . Concerning the fourth 
term in Eq. (Fl6l) . we show that 

JkiB jpDp q Bi m D mn diC n q J k iBi m B j n diD mn . 
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This property results from differentiating B iq — B im D mn C nq = with respect to z\ 
together with Eq. ([6]) and the antisymmetry of D. Using a circular permutation of 

k) these terms cancel with the JBBdD terms of Eq. ( |T5l) . 

After these steps, the Jacobi identity is rewritten as 

Sijk~\~Tijk~\~Uijk Bi m D mn Bl n BjpBf :: qOlDpq Bi m BjpBklD mn Dpq^ \C n q-\- (3(ijk) ■{l^J 

Next the strategy is to get rid of the dD terms. In order to do this, we insert C 
terms through a B coefficient in the first term of the previous equation, i.e., inserting 
Bj P = Bj a D a pCp p . From the identity 

CppdlDpqBkq = OlBk/3 — dlCfipDpqBkq — CppDpgdlBkq, 

which is obtained by differentiating Eq. §6§ with respect to Z\ (and using the 
antisymmetry of D and C), the first term in Eq. ( fTTl) is rewritten as 

Bi m D mn Bi n B jpBkqdiDpq Bi m D mn Bi n Bj a D a pDpqB] t qOiCfjp 

+ Bi m D mn Bi n Bj a D a p(diBki3 — CppDpqdiBkq). (18) 

From Eq. ([7j) we replace dC by VC. The first term becomes 

BimDmnBj a D a pBk q DqpBi n diCf3p+ 0(ijk)= B im D ' mn B j a D ' a pBk q D 'qpC n i^ 1 \C 'p p + O(ijk) ■ 

From the equation B im D mn C n i = Bu [see Eq. (jfjj)] we see that the previous term cancels 
with the second term in Eq. (1171) . still using a circular permutation of the indices k). 
In a similar way, the second term in Eq. f|T8|) vanishes by inserting Bj a D a pCj3 P = Bj p . 
Consequently, we have proved the Jacobi identity for the bracket (CQ) with Eqs. ([3)-®, 
i.e., 



3. Examples 

3.1. Example 1: Odd number of constraints 

First we describe a rather trivial example in order to illustrate the method. We consider 
the following Poisson matrix 



/ 





-Z3 


Z2 








\ 









-Zl 














Zl 



























-1 




V 











1 





/ 



which corresponds to the Poisson bracket 

dF x 8G dF dG + OF dG 
dz dz dw\ dw2 dw2 dw\ 
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We impose three constraints $fc(z, w) = Zk for k = 1,2,3. The associated operator Q 
is given by 

/ 1 \ 



Q 



\ 



10 
10 



The matrix C is a 3 x 3 antisymmetric, so it is not invertible. Its expression is 





/ 





-%3 


Z2 


\ 


c = 




Z3 





-zi 






V 


-z 2 


Zl 





/ 



The following matrix D satisfies Eq. 



D 



/ \ 
1 
0-10 



The solution of Eq. (j4]) is not unique, but all solutions lead to the same Dirac bracket 



J * 



which corresponds to the Dirac bracket 



UD: 












(° 














\ 
















































-1 




v° 








1 





/ 



{F, G}* 



OF dG OF dG 



dwi dw2 dw2 dwi 
3.2. Example 2: Vlasov-Poisson equation 

The second example concerns the Vlasov-Poisson reduction from the Vlasov-Maxwell 
equations. The field variables x(z) = (/(x, v), E(x), B(x)) where z = (x, v), and the 
equations of motion are 

/=-vV/-(B + vxB).a,/ ) 
E = V x B- J, 
B = -V x E, 

where J = J d 3 v v/. The Poisson bracket is given by 



{F, G} = I d G zF x -jG x , 
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where F x is the functional derivative of the functional F with respect to the field 
variables X- The Poisson matrix is given by (see Refs. [T ^ [T5 | [13]) 

/-[/,-] ~d v f \ 

J = -fdy 5(v)Vx , 
y -5(v)Vx j 

where the small bracket [•, •] is given by 

[f,g] = Vf ■ d v g -dvf-Vg + B- (d v f x d v g), 

where V (resp. d v ) is the partial derivative operator with respect to x (resp. v). 

In order to obtain the Vlasov- Poisson equations from the Vlasov-Maxwell equations 
we impose the following constraints: 

Q[f, E, B] (x) = (V x E, B - B (x)), (19) 

where B is a non-uniform background magnetic field. The operators Q is given by 



Q 



Vx 
1 



c 



(20) 



The operator C is given by 

(Vx) 2 
-(Vx) 2 

The operator C is not invertible; however, the Dirac procedure still applies with an 
appropriate choice for D given by 

A" 1 
-A" 1 



D 



(21) 



so that Eq. (Sj) is satisfied. We notice that Q(l - DC) ^ but - DC) = 0. This 

is due to the fact that V • B is a Casimir invariant of J. 

As a result, the Poisson operator of the Vlasov- Poisson equations is given by 



J * 



/ -[/,-] -VA^V-dv/ 0\ 

-VA^V ■ (fdy) 



V 















/ 



It leads to the expression of the Poisson bracket 

{F, G}* = J d 6 z f[F f - A- X V • F E ,G f - A- X V ■ G E ]. 

In the case of the constraints given by Eq. (fl9~|) . one of the constraints, V • B, is already 
a Casimir invariant of the Vlasov-Maxwell bracket, and hence a first-class constraint. 
This kind of redundancy is a source for the non-invertibility of the matrix C. In the 
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present example, we eliminate this redundancy by modifying the set of constraints (TT§j) . 
removing the divergence of B, i.e., we consider the following constraints : 

Q[f, E, B] (x) = (V x E, V(B - B (x))), 

where V is the projector on the solenoidal part given by V = 1 — VA _1 V-. Using similar 
calculations as above, we find 

Q=(° VX ° V 

yo o v J 

and C remains unchanged and is given by Eq. (|20|) . The matrix D given by Eq. (I2T!) 
satisfies Q)(l — DC) = 0. We have eliminated some redundancy at the origin of the 
non-invertibility of C. As we shall see below, this kind of redundancy is not an issue in 
the proposed approach to compute the Dirac bracket. 



4. Concluding remarks 

We consider an ensemble of M constraints $ n of which the first K ones are Casimir 
invariants of the original Poisson bracket (a particular family of first-class constraints). 
The matrix C is written as 

7 
C 



C 



where C is assumed to be invertible. Given that $^ for k = 1,...,K are Casimir 
invariants, JQ} = (0, JQ 2 ) where Q 2 is the matrix of the derivatives of . . . , $ n 

with respect to the phase space variables. It is straightforward to show that any matrix 
of the form 

-D\ 2 C- 1 

satisfies Eq. (|4]). Therefore in the definition of the constraints there is no need to worry 
about possible combinations giving rise to Casimir invariants of the original Poisson 
bracket. Notice that if some of the constraints are first-class but not Casimir invariants, 
this might lead to some inconsistencies in Eq. (Jl} like the ones in Refs. [T6| [6]. 

Another case of redundancy is when one (or several) constraints can be obtained 
from the other constraints. For instance, we assume that one constraint $i is dependent 
of the other constraints, i.e., $i = f(& 2 , . . . , <&m)- 111 this case, the usual Dirac's 
procedure cannot be carried out since the matrix C is not invertible. This matrix is 
given by 

' -b 1 " 
b C 



c 



where b = ({$ 2 , $i}, {$ 3 , $i}, . . . , {$m, and the coefficients of C are {$ n , $ m } 

for n, m > 2. In fact the vector (—1, d 2 , . . . , d^f) where d m = df/d<& m , belongs to the 
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kernel of C. We assume that C in invertible. Then the following matrix D satisfies 
Eq. ® 

' o -S 
d cr 1 



D 



where d = (d 2 , • • • , d M ). Here we have used the two properties, b^d = and Cd = b. 
As a consequence this kind of redundancy is properly handled using Eq. (Tj0) instead of 
the too stringent requirement that the matrix C is invertible. 

In summary the requirement for the existence of a Dirac-like bracket (0Q) is obtained 
by imposing Eq. (TjE]), i-e., that the constraints are Casimir invariants. It translates into 

Ker C C Ker jQ\ (22) 

which is a necessary and sufficient condition for the existence of a solution to Eq. (T4]). 
Indeed we denote r the rank of C and we rewrite it as 



C 



where O is an orthogonal matrix and C is invertible. A possible antisymmetric solution 
to Eq. (TjJ is given by 

°-o(\ 

given the condition ( 1221 ). 

If the specific choice of constraints is such that Eq. (1221) is satisfied, then the Dirac 
bracket (CQ) can be computed from Eq. (jlj), and it is a Poisson bracket everywhere in 
phase space, and not just on the surface defined by the constraints. Otherwise some 
obstructions are present, and one should modify the set of constraints so as to reduce 
the kernel of C. 
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